Real space investigation of structural changes at the metal-insulator transition in V0 2 
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Synchrotron X-ray total scattering studies of structural changes in rutile VO2 at the metal- 
insulator transition temperature of 340 K reveal that monoclinic and tetragonal phases of VO2 
coexist in equilibrium, as expected for a first-order phase transition. No evidence for any distinct 
intermediate phase is seen. Unbiased local structure studies of the changes in V-V distances through 
the phase transition, using reverse Monte Carlo methods, support the idea of phase coexistence and 
point to the high degree of correlation in the dimerized low-temperature structure. No evidence for 
short range V-V correlations that would be suggestive of local dimers is found in the metallic phase. 



PACS numbers: 71.30.+h 61.50.Ks 61.05.cf 

VO2 undergoes a transition from a metal to a non- 
magnetic insulator upon cooling below 340K.[T] Accom- 
panying this transition is a structural change from the 
high temperature tetragonal phase to a low temperature 
monoclinic phase, where pairing and tilting of vanadium 
ions result in chains with alternating long and short V-V 
distances along the c-axis (Fig. [l]).[2H!] The mechanism 
underpinning the insulating behaviour of the low tem- 
perature phase continues to attract widespread interest, 
and recent work suggests the simultaneous role of the 
V-V pairing and electron correlation to drive the sys- 
tem from a correlated metal to a slightly less correlated 
insulator. [5HI] Diverse methods have been used to under- 
stand this transition, including recent probes of dynamics 
and techniques sensitive to spatial inhomogeneities. [5HTT] 
The phenomenon as it occurs in nanostructures is also re- 
ceiving increasing attention. [T2HT4] 

Above 340 K, the rutile VO2 structure in space group 
P^/mnm has a single near- neighbor V-V distance of 




FIG. 1: (Color online) (a) The metallic high temperature 
tetragonal rutile form of VO2 with a single V-V distance and 
(b) the insulating low-temperature monoclinic form, showing 
dimerized chains of alternating short and long V-V distances 
along the c-axis. 



2.88 A.p~5j [16] Below 340 K, the monoclinic structure in 
space group P2\/c is characterized by V-V dimers along 
the c-axis displaying two distinct V-V distances: 2.65 A 
and 3.12 A.QMEE] The structures are displayed in Fig.[l] 
Interestingly, several recent papers address the formation 
of domain-like regions in monoclinic VO2 on approach- 
ing the transition temperature from below. [TJ HH [20] 
An unresolved question regarding the nature of these do- 
mains is whether they are associated with an intermedi- 
ate phase formed during the transition. A second ques- 
tion concerns the possible existence of dimers, albeit over 
short length scales, persisting in the metallic high tem- 
perature phase. 

Here we employ synchrotron X-ray total scattering in 
conjunction with real-space analysis — using least-squares 
as well as reverse Monte Carlo methods — to show that 
the transition from the low-temperature monoclinic to 
the high temperature tetragonal phase in rutile VO2 oc- 
curs in a first-order manner, with coexistence of the low 
and high- temperature phases at the transition tempera- 
ture. We employ pair distribution function (PDF) analy- 
sis, which has emerged in recent years as an indispensable 
tool for probing and understanding structural changes, 
and specifically, the range-dependence thereof. [2TH24] 
The advantages of the PDF method lie in its ability to 
probe structural changes that lack long range correla- 
tions, which are not captured by Bragg peaks in X-ray or 
neutron scattering. By analyzing PDF data with large- 
box reverse Monte Carlo (RMC) simulations, we are able 
to produce a supercell where the local environment is free 
from symmetry constraints of the high or low tempera- 
ture structures. [25H27] Atomic positions can therefore re- 
lax to best fit the experimental data, which includes both 
the Bragg and diffuse components. This approach pro- 
vides an opportunity to observe any continuum of atomic 
positions which might exist across the metal-to-insulator 
transition in VO2. 

Phase pure VO2 powders were prepared in evacuated 
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FIG. 2: (Color online) Least-squares fits to the PDF of bulk 
VO2 using (a) the low temperature monoclinic structure for T 
= 300 K (Rio, 11.6%) and (b) the high temperature tetragonal 
structure for T = 400 K (R w , 9.6%). Circles are experimental 
data and lines are fits to the average structure. The difference 
is displayed below each fit. (c) Changes in the R w fitting 
parameters, shown for cooling (o) and heating (x) through 
the transition temperature, show an abrupt transition with 
good fits at all temperatures outside of T = 340 K. 
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FIG. 3: (Color online) (a) Thermal parameters (JJij) ob- 
tained from PDF least-squares refinement as a function of 
temperature reveal abrupt changes at the transition temper- 
ature. Thermal ellipsoids (99%) for neighboring V polyhedra 
in the c direction are shown above and below the transition 
temperature, (b) V-V bond distances from the least-squares 
refinement (points) overlayed with a map of the V-V distances 
obtained from RMC modelling. 



and sealed silica ampoules from a stoichiometric mix- 
ture of V 2 5 (99.9%, Alfa Aesar) and V 2 3 , the lat- 
ter obtained by reduction of V2O5 in 5%-H2/N2 gas 
at 900° C. Synchrotron total scattering data on powders 
in Kapton tubes were collected in transmission mode 
at beamline 11-ID-B of the Advanced Photon Source 
at Argonne National Laboratory by using X-rays with 
energies near 90keV (wavelength A = 0.13702 A). Sam- 
ples were heated and cooled continuously from 250 K to 
400 K at a rate of 50 Kh -1 using an Oxford Cryosystems 
Cryostream 700. Scattering data were collected on an 
image plate system with a sample-to-detector distance of 
250 mm. Raw images were processed using fit2d.[28] 
PDFs were extracted as G(r) using pdfgetX2 with 
Qmax — 25 A -1 . [29] Least-squares profile refinements 
were carried out using PDFGUI.[30] RMC simulations 
were performed using rmcprofile[31] with a 10 x 10 x 16 
supercell, starting from atomic positions of the tetragonal 
structure. The simulations were constrained by the PDF 
(where G(r) above corresponds to D(r) as described by 
Keen [32 j) up to r = 25 A. RMC results shown in this 
work are averages of many simulations to ensure an un- 
biased representation of the fit to data. 

Least-squares refinements to all the PDFs at 2 K inter- 
vals were performed using both monoclinic and tetrago- 
nal structures [Fig. [2|a) and (b), respectively]. In the 
case of data fit to a monoclinic structure [Fig. J2ja)], we 
obtain a good fit from T = 250 K to 338 K. The high 



temperature tetragonal structure [Fig. J2^b)], results in a 
good fit to the data from T = 342 K to 400 K. The PDFs 
of the two structures are markedly different as seen in 
Fig.[2ja,b). There are no significant changes (other than 
lattice constants) in either end member until the tran- 
sition is reached. The sudden structural transition at T 
= 340 K is evident in plots of the goodness-of-fit param- 
eters R w to the high and low temperature structures as 
seen in Fig.[2|c). In addition to locating the transition 
itself, the evolution in both temperature regimes suggest 
pure, single-phase polymorphs of the VO2 end members. 
There is no broad hysteresis in the transition because the 
measurement allows thermal equilibrium to be reached at 
each point. In addition, the sample is a bulk powder of 
large crystallites, which implies the possibility of many 
concurrent nucleation events. 

Dimerization of the V cations leads to splitting of the 
2a Wyckoff site in the tetragonal VO2 structure. If 
a dimerized structure is fit using the high-temperature 
model, the structural changes must be accommodated 
by an increase in the V atomic displacement parameters 
Uij, which could convey information of the directional- 
ity of the displacement. Fig.[3|a) shows the least-squares 
refined values of the Uij parameters of the tetragonal 
structure over the full temperature range. At all tem- 
peratures, Uu = U22 by symmetry. Above the tran- 
sition, U33 is equal to Uu, while Z7 12 (the only non- 
diagonal parameter allowed by symmetry) is negligible. 
The high-temperature V position is therefore spherical, 
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and is shown on the right in Fig.[3ja). Below 340 K, 
dimerization along the c direction leads to a large in- 
crease in U33. Tilting of V-V vectors off the c axis (zero 
above the transition) produces an increase in Ui2- The 
corresponding growth and elongation of the V thermal 
ellipsoids are seen on the left in Fig.^a). 

The split V positions and corresponding dimerized V- 
V distances are the hallmark of cooling through the tran- 
sition, and we use the PDF to extract the distances di- 
rectly. If there were a gradual shift in V-V distances 
over the extent of a wide transition region, which other 
probes such as the electrical resistivity suggest, the lo- 
cal structure information in the PDF would reproduce 
these distances regardless of whether or not they are 
correlated over long ranges. Least-squares refinements 
of the V-V distances (using the appropriate monoclinic 
or tetragonal model for each structure) are displayed 
as points in Fig.|3jb) and show an abrupt change upon 
cooling through the transition. The distances at 336 K 
and 342 K are effectively unchanged from the low- and 
high-temperature structures at 300 K and 400 K respec- 
tively, so the structural transition occurs within this 
range (AT < 6K). The least-squares refinements, how- 
ever, are predicated on the choice of a monoclinic or 
tetragonal unit cell. 

Histograms of V-V distances from RMC supercells 
(not constrained by symmetry) are displayed as an inten- 
sity map in Fig.[3^b). They display excellent agreement 
with the least-squares refined values for all points except 
T = 340 K, implying that the single-phase least-squares 
refinements accurately reproduce the true distribution of 
V-V distances in the end members. However, at T = 
340 K the least-squares refined distances are of interme- 
diate length, while the RMC V-V histogram is unclear. 
Is there an intermediate phase with V-V distances that 
are distinct from the monoclinic and tetragonal struc- 
tures, or is there a mixture of two phases? Our subse- 
quent analysis of the linear combination of experimental 
PDFs, along with V positions from RMC simulations, 
can resolve this issue. 

In Fig. [4^ a), experimental PDFs are displayed for T = 
334 K, 340 K, and 342 K: below, during, and above the 
structural transition temperature. The T = 334 K and 
342 K PDFs can be refined to the end member struc- 
tures, so a linear combination of these two PDFs pro- 
duces a two-phase PDF containing three V-V distances: 
two monoclinic and one tetragonal. This linear combina- 
tion PDF is displayed as the points behind the T = 340 K 
data, and tracks the data up to 25 A. The difference curve 
between the combination PDF and the T = 340 K data is 
shown in Fig.[4^b), showing that the two PDFs are iden- 
tical within the limits of experimental noise. The agree- 
ment implies that, even on a local scale, the VO2 sample 
at 340 K comprises two structurally distinct regions, cor- 
responding directly to the tetragonal or the monoclinic 
structures. No additional contribution to the PDF is 
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FIG. 4: (Color online) (a) Experimental X-ray PDFs from 
V0 2 at T below (334 K), equal to (340 K), and above (342 K) 
the transition temperature. A linear combination of the 334 K 
and 342 K data is shown as points, along with a difference 
curve to the 340 K data. In (b) the full fitting range is shown- 
at the transition temperature the sample is comprised of two 
distinct end-member phases. 
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FIG. 5: (Color online) (a) The relative probabilities of V atom 
positions as obtained from RMC modeling are dimerized at 
four temperatures below and spherical at three temperatures 
above the transition. At T = 340 K, the probability distribu- 
tion is intermediate. In (b), a two-phase model is formed by 
averaging the experimental maps from T— 334 K and 342 K, 
both of which resemble the low- and high-temperature end 
members. This resulting map appears similar similar to the 
340 K data reproduced in (c). Scale bars are 0.2 A. The color 
bar indicates minimum and maximum probability. 



found. 

Corroborating evidence for the two-phase combination 
is seen in RMC simulations. We examine the distribu- 
tion of V positions around the ideal crystallographic site 
in real space by folding the positions of all V in the RMC 
supercell back into a single unit cell, which produces a 
cloud of 1600 V positions (from a 10 x 10 x 16 supercell) on 
each site. This cloud is viewed as a two-dimensional his- 
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togram showing the most probable positions of V cations 
in the a — c plane in Fig. [5] For each cloud with T > 
340 K, V distributions are spherical (as seen in the el- 
lipsoids in Fig.[3ja)) and centered on the ideal tetrago- 
nal positions. For T < 340 K, dramatic splitting of the 
V position is seen. These two spots correspond to the 
end member monoclinic positions. We see no evolution 
of V clouds versus temperature far from the transition, 
but the 340 K cloud possesses an intermediate shape, dis- 
tinct from the two regimes on either side. Just as was 
performed for the experimental PDFs themselves, we can 
linearly combine RMC results to produce models of a 1:1 
mixture of the two phases. In Fig.[5|b), an average of the 
clouds from 334 K and 342 K is displayed, as a model of 
coexistence of the end members. The experimental cloud 
from T = 340 K is reproduced in Fig.[5|c), and their simi- 
lar appearance indicates that this is a two-phase mixture, 
in agreement with the linear combinations of the PDF in 
Fig.g 

Qazilbash et a/. [9 have recently suggested the forma- 
tion of "metallic nanopuddles" at the transition temper- 
ature, and Cava and coworkers [7] have suggested the im- 
portance of metal-metal dimers in a series of substituted 
Yi- x M x 2 (M = Nb, Mo) samples. What has not been 
clear up to now is whether these implied a third, inter- 
mediate phase of VO2. Here, we have shown that at the 
transition temperature there are only two distinct phase 
populations: low temperature monoclinic VO2 together 
with the high temperature tetragonal phase. Above the 
transition, there is no evidence for the presence of V-V 
dimers. 
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